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Abstract 

I consider the next-to-leading high-energy cluster expansion of large-k double jet 
production in QCD, and I determine the corresponding one-loop quark and gluon 
impact factors for a self-consistent energy scale. The result shows that coherent 
angular ordering of emitted gluons holds for hard emission also, and singles out a 
scale which in essentially the largest virtuality in the process. Both remarks are 
relevant for the precise determination of the BFKL kernel at the next-to-leading 
level. 



PACS 12.38.Cy 



1 



Various attempts to understand the small- x HERA data on structure functions and 
final states have stimulated the analysis of high-energy QCD beyond the leading- 
logarithmic approximation ||. 

After several years of investigation, the high-energy vertices needed for the next-to- 
leading (NL) approximation have been computed and the irreducible part of the 
NL kernel of the BFKL equation || has been written out explicitly |9|-|TT|. It turns out 
that the NL terms provide (sizeable) corrections to the hard Pomeron intercept and to 
the singlet anomalous dimensions |H| . 

However, the definition of the NL kernel is not free of ambiguities which have pre- 
vented, so far, a complete quantitative analysis of the NL results. The most important 
ambiguity is due to the dependence of some NL features on the determination of the 
physical scale of the squared energy s, in the log s dependence of the cross sections. 

In a hard process like DIS, the scale of s is taken to be Q 2 , the virtuality of the photon 
or electroweak boson involved. Thus the structure functions are basically dependent on 
the Bjorken scaling variable x = Q 2 /s, with scaling violations induced by a s (Q 2 ). Similar 



considerations can be made for double DIS [12| or quarkonium production, where two 
hard scales are present. 

On the other hand, the high-energy cluster expansion needed for the definition of the 
NL kernel, has been mostly investigated in the case of parton-parton scattering, in which 
no physical hard scale is present. In order to define such scale (s), we have to deal with a 
partial cross section, by fixing the virtualities kj., k 2 of the off-shell gluon Green's function 
by high-energy factorization (Fig. la). Experimentally, this corresponds JL3] to the cross 
section associated with the production of two jets with given transverse momenta ki,k2 
and large relative rapidity, in the fragmentation regions of the incoming partons A, B. 

We can then define the scale So(ki,k 2 ) as the one occurring in the logarithmic be- 
haviour of the gluon Green's function G(ki, k 2 ; s), whose loop expansion reads 
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s (ki,k 2 ) ' """"" J 1 I 2"" u \~~°s ) ' " JV ^"°s 



G = 1 + I K Q log y~~ - r + const. ) + I }-K 2 ( log — ) + K NL log — + const. I 



(1) 



where we have introduced the leading and next-to- leading operator kernels Kg, Knl- It 
becomes then clear that the determination of K^l on the basis of an explicit two-loop 
calculation is dependent on the choice of sq: a change in the latter is a NL effect which 
induces a change in Knl, for the expansion ([I]) to be left unchanged. 

More precisely, I will assume in the following a k - factorized form for the differential 
cross section 

da AB f du ( s 



d[ki]d[k 2 ] 

where d\k] = rf 2 ( 1+e ^k/7r 1+<: in D = 4 + 2e dimensions, 

G U , = [1--(K Q + K NL )]- 1 (3) 

is the Mellin transform of the gluon Green's function [] and h A,B are (process dependent) 
impact factors, to be determined also, which may carry collinear singularities due to the 
initial (massless) partons A,B. 

In this note I will present arguments for a choice of Sq which is consistent with the 
collinear properties of the partonic process and I will give the corresponding one-loop 
determination of the impact factors. These results can then be used to yield a precise 
determination of Knl at two-loop level [fTOf] . 

It turns out that the basic physical issue is the so-called "coherent effect" in the space- 



like jets, noticed by the author long ago [ 14[ . This effect can be stated by saying that the 



leading log squared matrix element for gluon emission holds in a dynamically restricted 
phase space, in which angular ordering of the emitted gluons is required. The angular 
restriction determines, among other features, the energy scale also. 

Let me start noting that the one-loop Ag*{k. 2 ) — > Ag*(k 2 ) impact factor involves the 
known ]5| particle-particle reggeon vertex AA'g*{k.2) at virtual level, and the Ag*(k 2 ) — > 
A'g(q) squared matrix element (Fig. lb) at real emission level. The latter quantity can 
be extracted from the exact squared matrix elements of Ellis and Sexton [jl5|j , and has 
been investigated by various authors at matrix element level [|3],|8],[TB| but has not yet been 



^^An alternative formulation, not adopted here, could include an ^-independent NL factor in Eq.(3). 
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provided in explicit form in D dimensions, to the author's knowledge. I thus derive it 
directly for the quark case. 



By using the notation of Fig. 1 and the Sudakov parametrization 



k 2 

P1-P3 = K - z iPi - T\ + k i 

1 — Z1 IS 



Zi)S 

K = Zlrt-TT^YPt + K, (4) 

(1 - z 2 )s 

we want to compute the qAQB —> IaIbQ squared matrix element in the limit of large 
subenergy s 2 

s 2 = Zis > si = z i^) ^ ( q = ki + ks ) ) (5) 
zi{l - Z\) 

and fixed subenergy s\, so as to cover the parton A fragmentation region. Since s 2 is 
large, we can use k-factorization [|J to write 

where, by current conservation, we can use j| the Feynman gauge 

A ab = t a t b 1n(fa+ Ah\ _ t b t a l\{t>l- + 



^ 2p 3 q 2 Pl q 

1 

I*? 



[t a ,tb\-nr2T((ki - k 2 )^ x + (M+ A)g»x ~ (q + &i)a7/J • (7) 



By using physical polarizations with e • q = e ■ p 2 = and the high-energy kinematics 
we can write the colour decomposition 

A a ^en v \ = A ab = {t a , t b }A {+) + [t a , t b ]A { ~ ] , (8) 

with 

A ± = Ct(fa Pl -e- t\ fa A) + Ct\fap*-e+ l - £ jb fa) , 
4 1 _ 1 4 



By performing the straightforward gamma function algebra and the polarization sum 
in D = 4 + 2e dimensions, we obtain, for each colour structure, the result 



|A| 2 = (2s) 2 P q (e,z 1 ) (C 2 2p 3 q + C 2 (l-z 1 )2p 1 q + C l C 3 (2p 3 q(l-z 1 )+2p l q-z 1 k 2 2 )) 

(10) 

where the splitting function of q — > g type 

^-P gq = P q (e,z) = ±-(l + (l-z) 2 + ez 2 ) (11) 

appears to factor out in front of the transverse momentum dependence. Finally, by 
introducing phase space and coupling constant 

/ 1 \ 2 dz 1 d[k 1 ]d[k 2 ] g 2 s T(l-e) 
^ = H(^J * 2 *i(l-*i) ' tts = l4^' (12) 

and performing the colour algebra, we obtain 



do-V \ (2a s C F N e ) 2 P q (e, Zl ) 

= 7T- 



1 



N c a s (1 - z t )q. (q- ^k 2 ) 
7T q 2 (q-^ik 2 ) 2 



rffk^^]^;^ ' k?k 2 r(i- e 

7T q 2 (q- ^ik 2 ) 2 _ 

This result agrees with the one extracted from Ref.(15) in the high-energy limit. 

The cross section in Eq. (13) contains two colour factors, which have a simple interpre- 
tation, depending on the collinear singularities involved. The C\ term, with singularities 
at q 2 = ((q — -2ik 2 ) 2 = 0), comes from the Sudakov jet region, in which the emitted 
gluon is collinear to the incoming (outgoing) quark. 

On the other hand, the CfCa term is not really singular at either q 2 = or (q — 
z 1 k 2 ) 2 = 0, except for z\ = 0, which corresponds to the central region. It comes from the 
"coherent" region in which the gluon is emitted at angles which are large with respect 
to the qAqA' scattering angle, and is thus sensitive to the total qAQA' charge Ca ■ This is 
the region we are interested in, which is relevant for the energy scale, because it tells us 
how the leading matrix element, valid in the central region, is cut-off in the fragmentation 
region. 



Since the CfCa term in Eq. (13) acquires the 1/q 2 singularity for |q|(l — Z\) |ki|zi, 
we can write roughly, in the quark case, 



d\k.i]d\k. 2 ]dzi 



) = 7rhf\k 1 )hf\k 2 )^P q (e,z 1 )e(q(l-z 1 )-k l z 1 ) , (14) 

/ coherent 1 



hf(k) = 2a s C F N e /k\a s r- Xr "» 



7T 



The latter is the coherence effect prescription of Ref. (14), in which angular ordering 



7>rr.e.>^ ( 15 ) 

Z\ V — Z\ 



replaces the smooth behaviour of the differential cross-section in Eq. (13). The expression 
(14) is actually exact upon azimuthal averaging in ki of Eq. (13), at fixed |q|. 

Let me stress the point that in Eq. (14) the q —>■ g splitting function ~ Pq(e, z\) factors 
out in front of the leading matrix element for any value of Zi , not only for z\ <C 1 as 
assumed originally fl4j| . Therefore, Eqs (13) and (14) form the basis for the interpolation 
between Regge region [z,\ <C 1, any k 2 /q 2 ) and collinear region (k 2 /q 2 <C 1, any Zi), in 
small- x equations of CCFM type [0,|T7]. 



We are now in a position to start checking Eqs (1) and (2) and determining scale 
and impact factors. Let me start with a qualitative argument based on Eq. (14). By 
introducing the gluon (quark) rapidity y = \og(z\*/s / q)(Y\ = log(y/s / ki)) , the restriction 
(15) on half the phase space becomes 

< y < Y x - log f 9 9fcl = log(yi/Max(g, h)) . (16) 

In other words, the pure phase space q/y/s < Z\ < 1, (or < y < log(\/s/q)) is cut-off 
by coherence if q < ki, yielding q/y/s < Z\ < q/k%. 

This means that, by adding the remaining phase space related to the fragmentation 
region of B, the differential cross section is roughly 

da(1) \ -1.(0) fi . M.fOU,. A„„ 



d[ki]d[k 



2J/ coherent <T V Max(g, k x ) M&x{q,k 2 )J 



(17) 



Therefore, the physical scale for the energy is Max (q,ki)- Max (q,k 2 ), which can be 
roughly replaced by the expression 

So(k!,k 2 ) = Max(k2,k2) , (18) 

which has the same behaviour in the regions q <C k\ ~ k 2 , q — k 2 ^> k\ and q ~ k\ ^> k 2 . 

We have thus understood the mechanism by which the physical scale is provided by 
the largest virtuality, in qualitative agreement with what we know about deep inelastic 
scattering. 

The coherence argument above can be made more precise by using the exact squared 
matrix element in Eq. (13), in the phase space region connected with the fragmentation 
of A. In fact it is easy to check the representation 



dz 1 T-J7Z — r-z =71 — r log — — pr + h q (q, k X , 



q/v - s d^d^dzj CfCa r(l-c) q 2 V Max(g,A; 1 ) 

(19) 

where the appropriate scale has been subtracted, and 



/i,(q,ki) = / dzi 
Jo 



' (i-^i)q-(q-^ik 2 ) 1 , h 

(q-^iK 2 ) zi q 



(20) 



is a finite constant which, by Eq. (14), vanishes in the q — > /irmt, thus eliminating the 
g = singularity in front. 

The result (19), upon symmetrization, coincides with Eq. (17) for the logarithmic 
piece. Furthermore, the constant piece h q can be integrated over k x at fixed k 2 and 
interpreted as one-loop contribution to the quark A impact factor in Eq. (2): 



^ 0) (k 2 )A cal r(i- e ) 



/ ^Vq.ki)=^) e (-^-W)-^) • (21) 



The 1/e pole in Eq. (21) comes as no surprise, because it corresponds to the collinear 
divergence at k\ = due to the q — > g transition of the initial massless quark 

{2C F )- l % q {uj = 0) = jf 1 d* (p ff (e, *) - 1) = ~ + \e , (22) 



the - term being included in the leading kernel K . 

Eq. (19) yields, upon symmetrization, the scale Max (q,ki)- Max (9,^2) mentioned 
before. The latter can be translated into the scale Max (kf, k\) of Eq. (18) by the change 
in the constant piece 

5h q (q,ki) = log^G fc2fcl - log^e gfcl , (23) 

which vanishes at q = also and, upon ki integration, corresponds to the change in 
impact factor 

5h<p/hW = a s l ~i>\l). (24) 

Eqs (18), (19), (21) and (24) summarize the one-loop results for scale and quark 
impact factor at real emission level. The virtual corrections are known |JJ for both Cp 
and CfCa colour factors. The Sudakov Cp term turns out to cancel completely with the 
corresponding z\ - integrated real emission (u — moment) and is thus of little concern 
here. 

The CfCa term, on the other hand, yields || 

d ° ) = vr/ i f(k 1 )^ 0) (k 2 )2 W (fc 1 2 ) ( 5 2 ( 1 ^(q)vr 1 ^ • (25) 

/ virtual 



d[ki]d[k 



, s fll N F \ f nlU „ 85 5 Np , 



,o;(k 3 ) = -^k 



a Sl o £ r 2 (l + e) 



2e r(l + 2e) ' 



Here the logarithmic term provides the gluon trajectory renormalization, which regularizes 
the leading BFKL kernel, the second term in square brackets is a beta-function coefficient 
which yields the running coupling for the overall a 2 s factor, and half the last term adds 
up to Eq. (21) to yield the total one-loop correction to the quark impact factor, for the 
scale Max (q, ki)- Max (q, k 2 ): 



ft«(k)/<l(k)=a,(k 2 ^ 



4e 4/ 2 V 18 v 1 ; 9 N c 



(26) 



In the final result, Eq. (26), we have singled out the collinear singularity due to 
Eq. (22), so that its finite part is the last term in square brackets. 
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The above calculation can be repeated for initial gluons. If D = 4 the squared gg* — > 
gg matrix element can be extracted from the helicity vertices of Ref. (8) to yield 

da {1) \ n(2a s N c ) 2 N c a s (1 - z 1 )q(q - Zl k 2 ) + zfkj 

where 

P (z) ~ 1 + " 4 + ( 1 -^ 4 (28) 

PM) ~ 2z{l _ z) (28) 

is related to the g — > g splitting function. The Sudakov and coherent regions have in this 
case the same colour factor, due to the final gluons' identity. 

The main difference with the quark case comes from the non singular part of the 
splitting function, which yields a different constant piece in Eq. (19) and thus a different 
real emission contribution to the gluon impact factor. The latter has the expected collinear 
divergence related to 7 gg (^ = 0), and a finite term, whose precise determination requires 
generalizing Eq. (27) to 4 + 2e dimensions, and adding the gg* — > qq contribution. The 
scale determination remains, of course, unchanged. 

I thus conclude that the NL one-loop calculation determines the partonic impact 
factors with consistent collinear behaviour, depending on the choice of scale s — Max 
(kf,k|). With this choice the LL logarithmic term has no spurious q = singularities 
(as pure phase space would imply) and the NL constant piece has no q = singularity at 
all. Futhermore, no spurious kf = singularities are present (as, for instance, the scale 
s = kik 2 would imply). 

This means that the left-over ambiguity comes from a change of scale keeping the 
behaviour of Eq. (18) in all the regions k\ <^ k 2 ,k 2 <C hi, q <C k\ ~ k 2 as, for instance, 
so = Max (ki,q)- Max (k 2 ,q) does. In this case the change of scale can be reabsorbed 
in a change of impact factor without spoiling the collinear properties, as already done in 
Eq. (23) for the example just mentioned. 

Once a self-consistent scale is chosen and the corresponding impact factors are found at 
one-loop level, then the two-loop calculation will determine the NL kernel unambiguously. 
The corresponding results and physical implications will be reported elsewhere. 
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